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Abstract — A  novel  vapor  sensor  that  uses  polarization  interfer¬ 
ometry  in  a  form  birefringent  porous-silicon  film  is  introduced,  an¬ 
alyzed,  and  experimentally  characterized.  Simulations  and  anal¬ 
ysis  of  accuracy,  versatility,  stability,  and  control  of  dynamic  range 
of  the  device  are  provided.  The  simulation  accurately  predicts  the 
polarization  interference  signal,  which  is  used  to  estimate  the  ef¬ 
fective  refractive  indexes  characterizing  the  form  birefringence  of 
a  porous-silicon  film  with  0.001  accuracy.  The  device  was  tested 
for  the  detection  of  heptane  concentration  in  a  range  of  342-20  000 
ppm  (=2.0%). 

Index  Terms — Form  birefringence,  polarization  interferometry, 
porous  silicon,  vapor  sensor. 


I.  Introduction 

IMMENSE  EFFORTS  have  been  made  recently  in  devel¬ 
oping  optical  devices  made  of  porous  silicon  [1]  utilizing 
well-established  microfabrication  techniques.  Having  many  ad¬ 
vantages,  such  as  large  interaction  surface  area  for  sensing,  sim¬ 
plicity  and  repeatability  of  fabrication,  and  compatibility  with 
other  well-established  silicon  technologies,  porous  silicon  has 
attracted  much  attention  in  the  field  of  sensor  technology.  For 
hazardous  environment  applications,  such  as  sensing  flammable 
vapors  and  gases,  optical  sensors  attracted  appreciable  attention 
in  contrast  to  well-established  electronic  methods  based  on  ca¬ 
pacitance  [2]  and  conductivity  [3]  measurements.  There  exist 
several  different  methods  of  optical  sensing  based  on  optical 
reflectivity  [4],  [5],  birefringence  [6],  [7],  and  using  proper¬ 
ties  of  optical  waveguides  [8].  However,  their  practical  imple¬ 
mentations  are  somewhat  complicated,  and  the  characterization 
shows  some  discrepancies  between  calculations  and  measure¬ 
ments  [4]— [8]. 

In  this  manuscript,  we  report  a  novel  realization  of  a  vapor 
sensor  using  a  free-standing  porous-silicon  film.  The  method 
employs  a  form  birefringence  effect  in  porous-silicon  film  to 
construct  an  ultracompact  polarization  interferometer.  The  sen- 
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Fig.  1.  Schematic  of  experimental  setup  and  the  index  ellipsoid  of  a 
porous-silicon  film. 

sitivity  of  the  device  is  achieved  via  the  dramatically  enhanced 
interaction  area  of  the  porous  silicon  combined  with  high  sen¬ 
sitivity  phase  measurement  via  the  analysis  of  the  polarization 
interference  signal. 

II.  Experimental  Setup  and  Polarization 
Interferometry  With  Porous-Silicon  Film 

Porous-silicon  films  are  fabricated  in  a  heavily  doped  p-type 
<100>  silicon  wafer  by  anodic  etching  in  aqueous  high-fre¬ 
quency  (49  wt%)  solution  with  ethanol  in  the  volume  ratio  of 
3 :  1  at  room  temperature  [9].  Atomic-force  microscopy  was 
used  to  determine  the  average  pore  size  of  less  than  20  nm 
in  diameter,  whereas  the  porosity  of  82%  was  estimated  using 
the  gravimetric  method.  The  diameter  of  a  fabricated  porous- 
silicon  film  is  about  10-13  mm.  As  the  micro-pores  are  dis¬ 
tributed  uniformly  with  their  axis  perpendicular  to  the  surface 
of  silicon  wafer  [10],  the  resulting  artificial  form  birefringent 
porous-silicon  film  can  be  considered  to  be  a  positive  uniaxial 
anisotropic  material  due  to  the  in-plane  rotational  symmetry,  re¬ 
sulting  in  optic  axis  collinear  with  the  direction  of  the  pores. 
Consequently,  the  effective  dielectric  constant  for  the  optical 
field  polarized  in  the  plane  of  incidence  (i.e.,  transverse  mag¬ 
netic  (TM)  or  extraordinary  polarization)  will  depend  strongly 
on  the  angle  of  incidence  (see  Fig.  1),  whereas  the  effective  di¬ 
electric  constant  for  the  optical  field  polarized  perpendicular  to 
the  plane  of  incidence  (i.e.,  transverse  electric  (TE)  or  ordinary 
polarization)  will  remain  constant.  The  strong  polarization  dis¬ 
persion  in  a  porous-silicon  film  can  be  used  to  construct  an  ul¬ 
tracompact,  high  sensitivity  porous-silicon  polarization  interfer¬ 
ometer  with  a  special  emphasis  on  sensors  applications. 

A  schematic  diagram  of  our  experimental  setup  is  shown  in 
Fig.  1.  For  our  experiments,  we  use  an  external-cavity  tunable 
diode  laser  operating  at  a  wavelength  of  1 .52  //in  to  minimize 
the  absorption  losses  upon  propagation  through  silicon  (i.e.. 
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photon  energy  below  silicon  band  edge).  The  laser  beam  is 
first  introduced  into  a  polarizer  and  then  used  to  illuminate  the 
porous-silicon  sample.  The  transmitted  light  is  then  analyzed 
by  a  polarization  analyzer  and  introduced  into  a  photodetector 
for  measurements.  The  porous-silicon  sample  is  mounted  on  a 
precision  rotation  stage,  allowing  precise  control  of  the  surface 
normal  orientation.  Initially,  we  align  the  propagation  direction 
of  the  optical  beam  to  be  collinear  with  the  axis  of  the  pores, 
i.e.,  the  optic  axis,  setting  a  =  0.  For  sensing  experiments,  we 
insert  the  porous-silicon  sample  into  a  custom-made  glass  flow 
cell  connected  to  a  bubbler  for  the  exposure  to  solvent  vapors 
in  a  controlled  and  reproducible  manner. 

The  polarizer  and  analyzer  are  aligned  perpendicular  with  re¬ 
spect  to  each  other  and  at  45°  with  at-axis  in  the  x-y  plane. 
Under  simplifying  assumptions  of  minimal  absorption  and  sur¬ 
face  reflections  in  this  setup,  we  find  the  transmitted  intensity  at 
the  output  of  the  analyzer 

1(a)  =  ^  [ a 2  +  b2  —  2abcos(A(/)(a))]  (1) 

27r 

A <j>(a)  =  —  [«,,(«) /.,,((>:)  -  n0(a)L0(a)\  (2) 

where  a  and  b  are  the  amplitudes  of  transmitted  electric  fields 
for  TE  and  TM  polarizations,  respectively,  A <f>(a)  is  effective 
phase  difference  for  ordinary  and  extraordinary  polarized  com¬ 
ponents,  and  A,  rie(0)(o!),  and  Le^(a)  are  the  laser  wavelength 
in  vacuum,  angle-dependent  effective  refractive  indexes,  and  the 
geometric  path  lengths  in  a  film  for  extraordinary  (ordinary) 
component  for  given  angle,  respectively.  Note  that  n0(a)  =  n0 
(ordinary  index)  and  ne(a)  is  a  function  of  three  parameters,  a, 
na,  and  ne  (extraordinary  index)  in  an  index  ellipsoid.  Notice 
that  the  effective  phase  difference,  A<j)(a ),  does  not  depend  lin¬ 
early  on  the  angle  of  incidence  a ,  giving  rise  to  a  chirp-like  be¬ 
havior.  The  interference  signal  has  chirp-like  harmonic  depen¬ 
dence  on  ct,  which  is  important  for  accurate  analysis  and  better 
understanding  its  functional  dependence  using  numeric  simula¬ 
tion. 


III.  Experimental  Data  and  Analysis 

The  polarization  interference  signals  are  obtained  by  mea¬ 
suring  the  transmitted  optical  intensity  for  each  porous-silicon 
film.  Fig.  2(a)-(d)  shows  the  comparison  of  measured  and  sim¬ 
ulated  polarization  interference  signals  versus  the  angle  of  in¬ 
cidence  for  several  samples  of  porous-silicon  form  birefringent 
films  with  thicknesses  of 72, 100, 198,  and  228  fin i,  respectively. 

From  Fig.  2,  we  observe  a  chirp-like  harmonic  behavior  with 
increasing  chirp  as  the  thickness  of  the  form  birefringent  film 
increases  at  all  angles  of  incidence.  As  the  angle  of  incidence 
increases,  the  interval  of  the  oscillatory  signal  increases  (see 
Fig.  2).  This  peculiar  uneven  sinusoidal  variation  helps  the  data 
fit  process  to  extract  reliable  and  accurate  values  of  optical  in¬ 
dexes,  and  improves  the  fit  accuracy  of  all  extremums.  After 
obtaining  a  curve  fit,  any  slight  variation  of  the  indexes,  even  on 
the  order  of  ~  0.001 ,  deteriorates  the  accuracy  of  the  fit  signifi¬ 
cantly,  causing  some  of  the  peaks  to  shift  by  about  ~0.3°.  On  the 
other  hand,  if  both  index  values  ne  and  n0  are  varied  simultane¬ 
ously,  the  error  range  increases  [see  also  expression  for  A <j>(a) 
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Fig.  2.  Polarization  interference  signals  versus  angle  of  incidence  for 
porous-silicon  film  samples  with  thickness  of  (a)  72,  (b)  100,  (c)  198,  and 
(d)  228  ji  m.  Inset  figure  denotes  all  index  pair  coordinates  in  Table  I.  The 
simple  model  agrees  with  the  variations  of  all  the  index  pair  value  quite  well. 


TABLE  I 

Comparison  of  Optical  Index  Pairs  and  Error  Ranges  of  Four 
Different  Porous  Silicon  Samples,  of  Which  Porosity  is  About 
82%,  Estimated  by  Gravimetric  Method 


Porous  Si 

sample 

(rip,  n0)  lower 

(ne,  n0)  best 

(ne,  nD)  upper 

(±8ne,  ±8n0) 

72  pm 

(1.5330,  1.3788) 

(1.6394,  1.4533) 

(1.7393,  1.5283) 

(±0.10,  ±0.08) 

100  pm 

(1.5330,  1.3590) 

(1.6047,  1.4117) 

(1.6733,  1.4618) 

(±0.05,  ±0.05) 

198  pm 

(1.6733,  1.4635) 

(1.7066,  1.4877) 

(1.7393,  1.5117) 

(±0.033,  ±0.024) 

228  pm 

(1.5330,  1.3541) 

(1.5515,  1.3673) 

(1.5692,  1.3803) 

(±0.018,  ±0.013) 

in  (2)].  The  index  pair  values  of  best  fit  and  the  corresponding 
error  ranges  for  all  cases  of  transmitted  intensity  curves  in  Fig.  2 
are  summarized  for  comparison  in  Table  I. 

Note  that  the  error  range  decreases  as  the  thickness  increases. 
This  occurs  due  to  the  increase  of  phase  lag,  leading  to  more 
oscillation  periods  (zero-crossings),  thereby  improving  the  ac¬ 
curacy  of  the  fit.  It  is  interesting  that  the  coordinates  of  all 
index  pairs  in  Table  I  are  plotted  quite  closely  to  a  simple  model 
line  in  an  index  plane  (ne,  n0 ),  as  shown  in  the  inset  figure  of 
Fig.  2.  Although  detailed  discussion  on  this  model  is  beyond 
the  main  objective  of  this  letter,  this  simple  parametric  relation 
for  form-birefringent  dielectrics  can  be  derived  from  a  model  of 
a  simple  rectangular  cylinder  shape  pore  structure  with  a  para¬ 
metric  variable  related  to  porosity.  The  coincidence  between  the 
estimated  index  pair  values  and  a  simple  birefringent  model  is 
also  a  good  evidence  to  support  the  quality  and  reliability  of  the 
simulated  curve  fit  procedures  and  results. 
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Fig.  3.  Comparison  of  the  polarization  interference  signal  versus  angle  of 
incidence  for  a  porous-silicon  sample  of  thickness  198  /tm  in  air  and  in  vapor 
of  heptane  with  concentration  of  2.0%.  (a)  Enlarged  view  near  middle  peaks 
for  details,  (b)  Whole  range  view  for  the  comparison  of  pair  variations.  Note 
the  peak  shifts  are  not  linear. 

Next,  we  discuss  the  effects  in  porous-silicon  polarization  in¬ 
terferometer  that  lead  to  its  application  for  chemical  sensing. 
When  the  porous-silicon  film  sample  is  exposed  to  solvent  va¬ 
pors,  effective  indexes  nQ  and  ne  slightly  changed  due  to  the 
variation  of  surface  states  originated  from  the  difference  in  elec¬ 
tron  affinities  between  silicon  and  the  vapor  molecules  adsorbed 
on  the  surface  boundary  of  all  pores.  The  local  variation  of  index 
on  the  surface  boundary,  even  though  quite  small,  causes  shifts 
in  the  polarization  interference  signal  with  chirp-like  harmonic 
behavior,  and  thus,  it  provides  quite  large  change  in  the  output 
intensity  when  set  at  a  high  sensitivity  bias  angle  of  incidence. 
For  experiments  on  the  vapor  sensing  of  heptane  vapor  of  2.0% 
concentration,  the  change  of  the  polarization  interference  signal 
is  shown  in  Fig.  3  for  a  film  of  thickness  198  //m.  The  enlarged 
curve  view  is  shown  in  Fig.  3(a)  near  middle  peaks  for  details, 
and  the  whole  range  view  is  shown  in  Fig.  3(b)  for  the  compar¬ 
ison  of  curves  with  various  pair  variations.  Note  that  the  shift  of 
peaks  is  not  linear  in  a ,  although  it  is  not  so  clear  in  this  whole 
view. 

For  this  experiment,  we  used  same  sample  as  for  Fig.  2(c). 
The  pattern  for  “in  heptane”  was  obtained  by  a  small  adjust¬ 
ment  of  Sne,  as  shown  in  Fig.  3(a),  or  ( Sne ,  Sna),  as  shown  in 
Fig.  3(a),  in  the  simulations.  The  simulated  curves  in  entire  re¬ 
gion,  as  shown  in  Fig.  3(b),  enable  us  to  design  an  appropriate 
a  for  desired  sensitivity  by  revealing  all  nonlinear  shifts  in  a 
changed  pattern.  Note  that  if  the  pair  variation  exceeds  ( Sne , 
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Sn0)  =  (—0.02,  —0.0125),  the  fit  deviates  significantly  from 
the  measured  one,  while  three  curves  of  careful  pair  variations 
within  error  range  look  similar.  The  patterns  for  the  concentra¬ 
tion  of  heptane  from  342  to  20  000  ppm  (=2.0%)  are  fitted  well 
by  a  linear  variation  of  6ne.  Furthermore,  a  porous-silicon  film 
shows  quite  stable  behavior  for  repeated  experiments. 

IV.  Conclusion 

A  novel  vapor  sensor  using  a  freestanding  porous-silicon 
form  birefringent  film  is  realized  experimentally,  and  the 
characteristics  are  analyzed  using  computer  simulations.  It 
enables  us  not  only  to  design  an  appropriate  dynamic  range 
and  sensitivity  by  optimizing  the  angle  of  incidence  a,  but  also 
to  extract  accurate  index  values  of  birefringence  within  small 
error  of  just  ~0.001  for  index  variation.  Even  without  any 
careful  optimization,  the  concentration  of  heptane  vapor  from 
342  to  20000  ppm  (=2.0%)  can  be  detected  experimentally 
with  high  reproducibility.  A  freestanding  porous-silicon  form 
birefringent  film  provides  a  basis  for  the  construction  of 
ultracompact  polarization  interferometer  when  sandwiched 
between  a  polarizer-analyzer  pair,  enabling  fast  and  versatile 
sensing  mechanism  with  not  only  the  accuracy  from  the  high 
sensitivity  of  interferometry  but  also  the  convenience  from  the 
simplicity  of  implementation  for  the  detection  of  vapors  or 
gases,  especially  in  a  hazardous  environment. 
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